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Abstract.  The  relative  contribution  to  the  radiometric  force  of  the  area  and  perimeter  of  the  vane  is  studied 
experimentally  and  numerically.  Experimentally,  a  circle,  a  low-aspect  rectangle,  and  a  high-aspect  rectangle  were  all 
tested  on  a  thrust  measurement  stand  with  nano-Newton  resolution  placed  in  a  stagnant  gas.  The  computational  results 
were  obtained  through  2-D  simulations  using  the  direct  simulation  Monte  Carlo  method,  as  well  as  a  discrete  ordinate 
solution  of  the  ES  model  kinetic  equation.  Gas  pressure  was  varied  from  .006  to  6  Pa,  which  was  a  broad  enough  range 
to  observe  the  characteristic  peak  force  production  of  a  radiometer  in  the  transition  regime,  where  the  peak  occurs  at  Kn  ~ 
0.1.  It  was  found  that  the  area  of  a  radiometer  vane  is  responsible  for  a  significant  amount  of  the  total  force  production 
through  a  wide  range  of  its  operation.  It  is  only  at  the  highest  background  pressures,  well  after  force  production  has 
peaked,  that  the  vane  perimeter  appears  to  dominate  the  operation  of  the  radiometer. 
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INTRODUCTION 

The  objective  of  this  work  is  to  compare  experimentally  measured  forces  created  by  rarefied  gas  flows  on  heated 
plates  of  different  shapes  and  numerical  simulations  of  the  total  force  as  a  function  of  gas  pressure.  The  foundation 
for  this  work  was  laid  in  the  distant  past  by  William  Crookes,  1,2  who  invented  the  radiometer  in  the  1870’s.  It  was 
his  novel  apparatus,  and  his  pioneering  work,  which  helped  develop  and  validate  the  fledgling  branch  of  science 
now  known  as  gas  kinetics.  Simply  designed,  it  consisted  of  an  airtight  bulb  under  partial  vacuum  with  a  series  of 
half-blackened  vanes  mounted  on  a  spindle.  When  light  was  shone  on  the  device,  the  vanes  began  to  rotate  on  the 
spindle  with  the  half-blackened  sides  retreating  from  the  light  source.  Though  Crookes  initially  believed  he  was 
directly  measuring  photon  pressure,  a  series  of  investigations  by  other  prominent  scientists  of  the  era  proved 
something  entirely  different  was  responsible. 

Reynolds  first  suggested  that  the  force  was  due  to  surface  outgassing,  and  then  presented  a  more  rigorous 
explanation  based  on  kinetic  theory 3.  According  to  the  latter  theory,  the  gas  in  the  partially  evacuated  bulb  is  the 
main  driving  force  responsible  for  the  rotation  of  the  vanes.  Further  work  was  done  by  Schuster4,  which  turned  up 
the  first  experimental  evidence  of  gas  forces  being  the  dominant  cause  of  the  radiometric  effect.  The  kinetic  theory 
explanation  is  in  fact  a  free  molecule  approximation  of  the  radiometric  effect:  after  the  light  heats  the  blackened 
side  of  the  vane,  the  molecules  reflect  from  the  hot  side  with  higher  average  velocities.  This  leads  to  a  larger 
momentum  imparted  to  the  hot  side,  and  results  in  the  motion  of  the  vanes  with  the  hot  side  trailing. 

Reynolds  also  proposed  that  increasing  the  pressure  in  the  bulb  would  dramatically  affect  the  device’s  operation. 
He  theorized  that  the  onset  of  intermolecular  collisions  in  the  gas  would  work  to  equalize  the  pressure  on  both  sides 
of  the  vane,  and  under  sufficiently  high  pressures  (near-continuum)  the  force  production  would  drop  rapidly.  5  With 
further  exploration,  Maxwell6  showed  that  an  unbalanced  force  exists  near  the  edge  of  the  heated  side  of  the  vane, 
where  the  heat  flow  in  the  gas  is  non-uniform.  Though  a  unique  result  for  the  time,  it  took  another  fifty  years  before 
Einstein  presented  a  simple  theory7  that  related  the  force  on  the  vanes  to  their  perimeter.  This  edge  dependence  of 
the  vane  force  has  found  partial  confirmation  in  experimental  work8,  where  the  force  was  found  to  depend  on 
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perimeter.  Although  the  dependence  was  not  shown  to  behave  exactly  as  Einstein  had  predicted,  the  edge  theory  of 
the  radiometer  has  become  widely  accepted. 

It  was  with  this  understanding  that  interest  in  the  radiometer  problem  declined  steadily  after  1928,  mostly 
because  the  issue  of  force  production  was  considered  closed,  and  no  direct  application  for  radiometric  forces  had 
been  identified  at  that  time.  The  phenomenon  that  drives  Crookes'  radiometer  has  been  summarized  by  Draper9  who 
in  effect  described  our  present  understanding  of  it:  a  temperature  gradient  along  a  surface  creates  tangential  stresses 
in  the  gas,  and  as  a  result  the  gas  moves  over  the  surface  from  the  cold  to  the  hot  side.  For  the  gas  to  be  accelerated 
in  this  manner,  it  is  required  that  an  equal  and  opposite  force  propel  the  plate  in  the  opposite  direction.  Under  this 
explanation,  the  main  contribution  of  the  force  to  the  rotation  of  the  vanes  comes  from  the  an  area  very  near  the 
edge(a  zone  with  a  width  of  a  mean  free  path  X,  according  to  Einstein). 

Lack  of  an  obvious  application  may  explain  the  miniscule  interest  in  the  problem  in  the  second  half  of  the  last 
century,  when  only  a  few  research  papers  on  the  subject  were  published  10,  with  other  publications  being  historic 
analyses  and  overviews.  The  situation  has  changed  over  the  last  decade,  as  the  radiometric  phenomenon  has  been 
found  to  have  possible  uses  in  a  number  of  different  micro-  and  macro-scale  devices.  The  forces  produced  by 
temperature  gradients  in  a  gas  have  been  shown  11  to  be  applicable  to  modern  microactuator,  where  the  direct 
simulation  Monte  Carlo  (DSMC)  method  was  used  to  model  forces  on  vanes  mounted  on  an  armature.  This  method, 
along  with  experimental  measurements,  has  been  employed  12  to  study  a  concept  of  an  opto -microengine  that  uses 
radiometric  forces.  Subsequently,  a  series  of  papers  by  Passian  with  co-workers  was  published,  13, 14  where 
radiometric  phenomena  were  studied  experimentally  and  analytically,  mostly  with  application  to  microcantilevers. 
The  use  of  radiometric  forces  as  an  approach  to  study  gas-surface  translational  energy  accommodation  has  also  been 
mentioned  by  Passian  et  al  in  a  third  paper,  15  and  a  new  concept  of  a  high-altitude  aircraft  supported  by  microwave 
energy  that  uses  radiometric  effects  was  put  forward  by  Benfordl6. 

The  new  studies  of  the  old  radiometric  phenomena  have  been  supported  on  the  one  hand  by  modern  technologies 
that  allow  more  accurate  measurements,  and  on  the  other  hand  by  state-of-the-art  numerical  methods  that  rely 
heavily  on  parallel  computing.  These  two  factors,  along  with  the  revived  interest  in  the  application  of  radiometric 
phenomena,  have  prompted  the  authors  to  revisit  the  contribution  of  the  “collisionless”  (area)  versus  the 
“collisionar”  (edge)  forces  to  the  total  radiometric  force. 

EXPERIMENTAL  SETUP 

To  study  the  role  of  area  and  edge  forces  on  the  plate,  three  radiometer  vane  geometries  were  used.  Each  of  these 
consisted  of  a  Teflon  insulator  sandwiched  between  two  aluminum  plates.  A  resistive  heater  was  located  between 
one  of  the  plates  and  the  Teflon  insulator,  and  the  temperature  of  the  hot  side  of  the  device  was  maintained  by 
varying  the  power  input  to  the  heater.  Each  of  the  plates  and  insulator  are  0.32cm  thick,  and  when  assembled  yield  a 
total  device  thickness  of  0.95cm.  Figure  1  shows  a  CAD  model  of  one  of  the  devices  used  for  both  experimental 
and  computational  comparisons.  The  first  device  was  a  rectangle  3.81  x  12.7cm,  the  second  device  was  a  rectangle 
with  exactly  double  the  area  having  dimensions  of  7.62  x  12.7cm,  and  the  third  was  a  circle  with  a  diameter  of 
1 1.13cm  with  an  area  equivalent  to  that  of  the  large  rectangle. 

Each  of  these  devices  was  individually  mounted  on  a 
modified  nano-Newton  Thrust  Stand  (nNTS)17  located  inside 
of  a  0.4m  diameter  vacuum  chamber.  When  calibrated  using 
a  set  of  electrostatic  combs18,  the  nNTS  provides  very 
accurate  and  repeatable  data  with  typical  force  resolution  of 
approximately  0.1  uN  and  statistical  scatter  around  1%.  The 
experimental  error  based  on  standard  deviation  ranges  from  a 
few  percent  at  the  lowest  pressures  to  approximately  1% 
through  most  of  the  curve.  However,  due  to  the  normalization 
by  experimental  temperature  measurements,  the  variation  in 
background  pressure,  and  the  small  uncertainty  of  the 
calibration  method,  the  total  absolute  experimental 
uncertainty  is  ~4%.  Day-to-day  variation  of  multiple  data  sets 
has  been  observed  to  be  ~2%. 

The  experimental  data  was  obtained  for  each  device  by 
evacuating  the  chamber  to  a  base  pressure  below  10‘3  Pa.  A 


Teflon  Hot  side 


Cold  side  Heater 

FIGURE  1.  CAD  drawing  of  a  radiometric  device 
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constant  voltage  was  applied  to  the  heater,  and  this  resulted  in  the  main  surfaces  reaching  temperatures  of 
approximately  415K  (hot)  and  365K  (cold).  The  background  pressure  inside  the  chamber  was  varied  by  flooding 
the  chamber  with  a  specific  gas  over  a  range  of  pressures  from  0.1  Pa  to  1.2  Pa.  Argon,  helium,  and  nitrogen  were 
all  utilized  with  the  additional  use  of  xenon. 

The  experimentally  measured  force  was  normalized  by  the  temperature  difference  between  the  hot  and  cold 
plates,  AT,  for  the  purpose  of  comparing  results  for  different  pressures  and  gases.  Verification  of  the  validity  of  the 
experimental  normalization  method  is  demonstrated  in  a  prior  paper19,  where  exceptional  linearity  with  temperature 
difference  is  observed. 


COMPUTATIONAL  SETUP 

Two  kinetic  approaches  were  used  in  the  work,  the  direct  simulation  Monte  Carlo  (DSMC)  method  and  a  finite 
difference  solution  of  the  ellipsoidal  statistical  (ES)  model  kinetic  equation.  The  DSMC  computational  tool, 
SMILE20,  was  used  in  all  DSMC  computations  presented  in  this  work.  The  variable  hard  sphere  model21  was  used 
for  the  intermolecular  potential,  and  the  Larsen-Borgnakke  model  with  variable  rotational  relaxation  number  used 
for  translational-internal  energy  transfer.  The  gas -surface  model  was  assumed  to  be  diffuse  with  full  energy  and 
momentum  accommodation.  A  two-dimensional  module  of  SMILE  was  used  in  this  work  because  three-dimensional 
modeling  is  prohibitively  expensive  for  most  pressures  under  consideration.  The  computations  were  conducted  for 
two  plates  with  dimensions  0.95x3.8 lcm  and  0.95x7. 62cm  immersed  in  an  initially  uniform  stagnant  gas  in  a 
chamber.  The  third  (Z  axis)  dimension  of  12.7cm  was  assumed  when  calculating  forces  consistent  with  the  actual 
size  of  the  plate  in  the  present  experimental  study.  The  main  surfaces  of  the  plate  were  set  to  41  OK  (cold)  and  450K 
(hot),  and  the  sidewalls  were  430K.  The  chamber  wall  temperature  was  assumed  to  be  constant  at  300K.  The 
computations  were  performed  for  chamber  pressures  ranging  0.006Pa  to  6Pa  and  chamber  sizes  from  0.45cm  to 
1.8m  using  three  gases:  helium,  nitrogen,  and  argon.  For  the  solution  of  the  ES  model  equation,  argon  and  helium 
were  considered  with  the  same  flow  conditions  as  in  DSMC.  Three-blocks  non-uniform  rectangular  grid  was  used  in 
the  solution  of  the  ES  equation  with  the  total  of  13,359  spatial  nodes.  The  Gauss-Hermite  quadrature  of  order  16  and 
three-eighths  rule  were  used  for  discretization  of  velocity  magnitude  and  angle,  respectively. 

RESULTS 


FIGURE  2.  Flowfield  around  a  radiometric  device  FIGURE  3.  Pressure  on  main  surfaces  of  the  radiometer 

The  flow  field  structure  typical  for  the  transitional  regime  around  a  plate  is  shown  in  Fig.  3,  where  the 
translational  temperature  and  streamlines  are  shown  for  helium  at  a  chamber  pressure  of  2Pa.  The  computation 
domain  of  0.45x0.45m  is  bounded  by  300K  diffuse  walls.  The  hot  side  is  on  the  left  hand  side  of  the  plate.  There  are 
four  vortices  created  by  the  temperature  gradients,  two  at  each  side  of  the  plate.  The  vortices  at  the  hot  side  of  the 
plate  are  noticeably  stronger  than  at  the  cold  side.  For  the  case  under  consideration,  the  maximum  flow  speed  in  the 
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hot  surface  vortices  is  about  5m/s,  and  in  the  cold  surface  vortices  it  is  about  2m/s.  Note  that  the  vortex  structure  is 
similar  to  that  calculated  by  Ota12  but  qualitatively  differs  from  the  flow  structure  given  by  Kennard22 

To  examine  the  contribution  of  the  edge  effects  to  the  net  force,  and  analyze  the  relative  importance  of  the  forces 
on  the  edges  of  plate  versus  the  area  forces,  the  computations  were  performed  in  helium,  for  a  geometric 
configuration  of  a  solid  3.9cm  plate.  Consider  first  the  surface  pressure  distribution  over  the  cold  and  hot  sides  of  the 
solid  plate,  shown  in  Fig.  7  for  a  chamber  pressure  of  2Pa. 

The  pressure  has  a  minimum  near  the  center  of  the  plate  both  for  the  hot  and  cold  sides.  The  values  near  the 
edges  are  visibly  larger  than  that  near  the  center.  The  net  force,  however,  is  produced  by  the  difference  of  pressure 
forces  on  the  hot  and  cold  sides  of  the  plate,  and  the  contribution  of  area  and  perimeter  forces  may  be  compared 


analyzing  the  pressure  difference  near  the  center  of  the  plate  and  closer  to  the  edge. 
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FIGURE  4.  DSMC  (top)  and  ES  (bottom)  solutions  for 
He  at  2Pa 
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FIGURE  5.  Pressure  difference  along  the  vane  for  Ar 


The  statistical  scatter  inherent  in  the  DSMC  method  significantly  complicates  such  an  analysis,  and  the 
deterministic  solution  of  the  ES  model  kinetic  equation  was  therefore  used  for  this  purpose.  Although  it  uses  an 
approximation  to  the  collision  operator  of  the  Boltzmann  equation,  the  ES  equation  nevertheless  is  capable  of 
capturing  the  main  features  of  the  radiometric  flows.  Comparison  of  the  DSMC  and  ES  solutions  is  presented  in 
Fig.  4,  where  the  temperature  fields  and  streamlines  are  shown.  The  numerical  results  obtained  by  these  two 
approaches  are  similar,  and  the  difference  in  temperature  is  less  than  1  percent. 

The  pressure  difference  distribution  alone  the  vane  is  shown  in  Fig.  5  for  pure  argon  at  four  pressures.  Note  that 
the  pressure  of  1.219Pa  corresponds  to  the  maximum  radiometric  force.  Only  half  of  the  plate  is  shown  here,  and  the 
coordinate  Y=0  corresponds  to  the  center  of  the  vane.  For  the  lowest  pressure  shown  in  this  figure,  the  pressure 
difference  near  the  center  is  close  to  that  near  the  edge,  which  indicates  that  the  area  force  is  the  main  contributor  to 
the  radiometric  force.  For  the  higher  pressure,  the  opposite  is  true:  the  perimeter,  or  edge,  force  is  the  determining 
factor.  Both  area  and  perimeter  forces  are  contributing  for  gas  pressures  corresponding  to  the  maximum  radiometric 
forces,  with  the  area  force  appearing  to  be  somewhat  larger. 

Experimental  results  for  nitrogen  in  the  small  chamber  are  shown  in  Fig.  6.  It  is  readily  observed  in  the  low 
pressure  region  that  the  circular  plate  and  large  rectangular  plate  collapse  onto  each  other  and  that  their  values  are 
very  close  to  twice  that  of  the  smaller  plate.  This  is  expected,  as  free  molecular  theory  shows  the  force  to  be  entirely 
dependent  on  area.  As  the  flow  transitions  from  the  collisionless  regime,  the  picture  becomes  distinctly  more 
clouded.  While  it  is  readily  observed  that  the  plates  with  larger  area  produce  more  force  at  their  respective  peaks, 
the  force/area  ratio  does  not  hold. 

When  comparing  the  peaks  of  the  large  and  small  rectangles,  it  is  found  that  the  small  plate  creates  65%  of  the 
force  of  the  larger  one.  Interestingly,  this  is  almost  precisely  half  way  between  the  area  ratio  (50%)  and  the 
perimeter  ratio  (81%).  Were  that  the  only  data  presented,  it  would  be  tempting  to  presume  that  the  perimeter  force 
was  slowly  beginning  to  dominate.  This  does  not  entirely  appear  to  be  the  case  however,  as  the  peak  force  created 
by  the  circle  is  decidedly  larger  by  -7%.  It  seems  quite  contradictory  to  accepted  theory  that  a  device  with  14%  less 
perimeter  would  create  a  larger  force  without  an  equivalent  increase  in  heated  surface.  This  increase  in  force  is 
possibly  explained  by  more  efficient  pressure  redistribution  in  front  of  the  plate  due  to  the  axial  symmetry  of  the 
circular  plate  vortices. 
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The  numerical  predictions  obtained  for  two  different  plate  sizes  are  shown  in  Fig.  7  and  are  qualitatively  similar 
to  the  above  measurements.  In  the  free  molecular  regime  the  larger  plate  produces  a  proportionally  larger  force. 
The  ratio  of  the  small  plate  force  to  the  large  plate  force  increases  with  pressure  and  reaches  ~0.7  in  the  range  of 
pressures  where  the  force  is  maximum.  This  furthers  the  case  that  the  area  is  still  a  large  contributor  to  the  peak 
force,  where  it  must  be  noted  that  the  perimeter  ratio  of  the  2-D  case  is  unity. 
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FIGURE  6.  Normalized  force  for  different  device  geometries  in  the  0.4m  chamber 


FIGURE  7.  Numerical  results  for  two  different  plates  FIGURE  8.  Comparison  of  experimental  and 

computational  results 

Comparison  of  the  DSMC  results  and  experimental  measurements  of  the  net  force  on  the  plate  in  nitrogen  is  shown 
in  Fig.  8  for  different  pressures.  The  force  is  normalized  by  the  temperature  difference  between  the  cold  and  hot 
surfaces,  since  the  actual  values  of  the  temperature  somewhat  varied  in  the  experiment.  The  comparison  shows  that 
both  measurements  and  DSMC  predict  pressure  maximum  around  IPa.  The  numerical  results  noticeably  overpredict 
the  measurements,  with  the  difference  attributed  primarily  to  the  three  dimensionality  of  the  experiment  and  an 
incomplete  energy  accommodation  that  was  not  accounted  for  in  DSMC. 


CONCLUSION 


A  study  of  the  radiometric  forces  on  heated  plates  has  been  conducted  both  experimentally  and  computationally. 
The  experiments  were  carried  out  at  USC  in  a  0.4m  vacuum  chamber  up  to  a  maximum  pressure  of  1.2  Pa  for 
various  carrier  gases.  The  computations  were  performed  with  the  DSMC  method  and  a  solution  of  the  ES  model 
kinetic  equation  for  a  2-D  gas  flow  over  a  comparable  range  of  pressures.  It  is  shown  that  the  radiometric  devices 
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provide  maximum  force  at  a  Knudsen  number  approximating  0.1.  Qualitatively,  the  experimental  data  and 
computational  results  are  similar.  A  lack  of  experimental  data  on  gas -surface  accommodation  and  flow  three- 
dimensionality  yields  approximately  a  25%  difference  in  the  magnitude  of  the  measured  and  computed  forces. 
Comparison  of  three  geometric  configurations  has  shown,  both  experimentally  and  numerically,  that  the  effect  of  the 
area  is  significant  at  pressures  up  to  where  the  force  is  maximum. 
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